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Abstract 

The spin ordering and spin-charge coupling in LuFe204 were investigated on the basis of den- 
sity functional calculations and Monte Carlo simulations. The 2:1 ferrimagnetism arises from the 
strong antiferromagnetic intra-sheet Fe 3+ -Fe 3+ and Fe 3+ -Fe 2+ as well as some substantial antifer- 
romagnetic Fe 2+ -Fe 3+ inter-sheet spin exchange interactions. The giant magnetocapacitance at 
room temperature and the enhanced electric polarization at 240 K of LuFe204 are explained by 
the strong spin-charge coupling. 

PACS numbers: 75.80.+q,71.20.-b,77.80.-e,64.60.De 
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Recently, multiferroics |l|, |2j, y, |4j, |5|, la, 13, la, |9|, |10| have attracted much attention be- 



cause of their potential applications in novel magnetoelectric and magneto-optical devices. 
Among the newly discovered multiferroics, LuFe204 is particularly interesting due to its 
arge ferroelectric (FE) polarization |3j and giant magnetocapacitance at room temperature 
4J. In the high-temperature crystal structure of LuFe 2 04 with space group R3m, layers of 
composition Fe2C>4 alternate with layers of Lu 3+ ions, such that there are three Fe2C>4 layers 
per unit cell. Each Fe2C>4 layer is made up of two triangular sheets (hereafter, T-sheets) of 
corner-sharing Fe0 5 trigonal bipyramids (Fig. [1]). Below 320 K (T C o) LuFe 2 4 undergoes a 
three-dimensional (3D) charge ordering (CO) (2Fe 2 ' 5+ =>• Fe 2+ + Fe 3+ ) with the v^3 x \/3 
superstructure in each T-sheet; in each Fe204 layer, one T-sheet has the honeycomb net- 
work of Fe 2+ ions with a Fe 3+ ion at the center of each Fe 2+ hexagon (hereafter, the type 
A T-sheet), while the other T-sheet has an opposite arrangement of the Fe 2+ and Fe 3+ ions 
(hereafter the type B T-sheet). 

LuFe204, with the novel CO-driven "electronic ferroelectricity" , |3| presents several fun- 
damental questions. First, LuFe204 shows strong Ising behavior with the easy axis along c 
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12] . The spin anisotropy of the non-CO state is understandable because the spi n down 



electron of the Fe 2,5+ ion partially occupies the degenerate (d x 2_ y 2,d xy ) orbitals [5|, |l3j. How- 
ever, the Ising behavior below Tco is puzzling because the insulating \/3 x CO breaks 
the 3-fold rotational symmetry hence lifting the degeneracy of the (d 7 .2_ y 2,d xv )^ orbitals 



Second, LuFe204 undergoes a ferrimagnetic spin ordering below 240 K (Tn) [111, 
A number of experimental studies found this spin ordering to be two-dimensional (2D) in 
17| . In contrast, a recent neutron diffraction study observed a finite spin 
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correlation along c and suggested a 3D spin structure without considering CO [16|. The 



Mossbauer [14| and neutron diffraction 



15| studies led to a detailed ferrimagnetic structure 



of LuFe204, in which the majority spin lattice consists of all Fe 2+ ions plus one-third of the 
total Fe 3+ ions while the minority spin sublattice consists of the remaining Fe 3+ ions. This 
2:1 ferrimagnetic order was suggested to originate from weak ferromagnetic (FM) interac- 
tions between the next-nearest neighbor (NNN) Fe sites in the triangular antiferromagnetic 
(AFM) Ising lattice However, using the spin exchange parameters estimated from the 
energy parameters of LaFe03, Naka et al. [18| predicted quite a different spin structure that 
includes some Fe sites without unique spin direction. Therefore, the detailed ferrimagnetic 
structure and its origin remain unclear. Third, LuFe204 exhibits a giant magnetodielec- 
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trie response at room temperature [4J, and a room-temperature dynamic magnetoelectric 
coupling was also reported [l9]. Furthermore, the FE polarization of LuFe2C>4 was found 
to increase around T N jsj] . These observations suggest the occurrence of coupling between 
the CO and magnetism. The understanding of the spin-charge coupling is crucial for future 
magnetodielectric applications of LuFe204. 

In this Letter, we explore these isuues on the basis of first principles density functional 
calculations for the first time. A large spin anisotropy is found along the c direction due 
mainly to the Fe 2+ ions of the B-sheet, the spin ground state of the \/3 x \/3 CO state has 
the 2:1 ferrimagnetic spin arrangement proposed by Siratori et al. Jjj], and there occurs 
strong spin-charge coupling in LuFe 2 4 . 

Our density functional theory calculations employed the frozen-core projector aug- 
mented wave method [2[J encoded in the Vienna ab initio simulation package 21], and 
the generalized-gradient approximation (GGA) 22|. To properly describe the strong elec- 
tron correlation in the 3d transition-metal oxide, the GGA plus on-site repulsion U method 
(GGA+U) [23( was employed with the effective U value (U e ff = U — J with J = 0) of 
4.61 eV 5] . It is known experimentally , [l3] that the interlayer magnetic interactions 

in LuFe204 are weak, which is understandable due to its layered structure. In this work, 
therefore, we focus on the 2D spin ordering within a sing le Fe 2 4 layer. For the ^3 x 
CO state of LuFe204, the FE ordering of the Fe204 layers will be assumed. 

We first examine the magnetic anisotropy of the Fe ions by performing GGA+U cal- 
culations, with spin-orbit coupling (SOC) included, for the FM state of LuFe204 with the 
y/3 x v3 CO. As shown in Fig. QJa), there are two kinds of Fe 2+ ions and two kinds of 
Fe 3+ ions in the v3 x \/3 CO state. We label the Fe 2+ and Fe 3+ ions of the type A T-sheet 
as 2A and 3A, respectively, and those of the type B T-sheet as 2B and 3B, respectively. 
In our GGA+U+SOC calculations with spins pointing along several different directions, all 



Fe 2+ and Fe 3+ spins are kept in the same direction. Our ca 



axis is along the c direction, as experimentally observed 11 



culations show that the easy 



121 ] : the ||c-spin orientation is 



more stable than the _Lc-spin orientation by 1.5 meV per formula unit (FU). The orbital 
moments of 2A, 2B, 3A, and 3B for the ||c-spin orientation are 0.101, 0.156, 0.031 and 0.035, 
respectively, which are greater than those for the _Lc-spin orientation by 0.019, 0.062, 0.015, 
and 0.018 \Ib, respectively. As expected, the Fe 3+ (d 5 ) ions have a very small anisotropy, 
However, two kinds of the Fe 2+ ions also have different degree of spin anisotropy. The spin 
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down electron of the 2B Fe 2+ ion occupies the (d x 2_ y 2,d xy ) manifold j^], therefore the 2B Fe 2+ 
ion has the largest spin anisotropy along c. Our calculations indicate a non-negligible or- 
bital contribution to the total magnetization, in agreement with the X-ray magnetic circular 
dichroism result [l^ . 

To determine the magnetic ground state of LuFe 2 4 in the \/3 x y/3 CO state, we extract 
its spin exchange parameters by mapping the energy differences between ordered spin states 
obtained from GGA+U calculations onto the corresponding energy differences obtained from 
the Ising Hamiltonian 241 ]: 

H = JijSi Z Sj Z , (1) 

where the energy is expressed with respect to the spin disorder (paramagnetic) state, Jjj is 
the spin exchange parameter between the spin sites i and j, and S,- lz is the spin component 
along the c direction (\S Z \ = 2 and 2.5 for Fe 2+ and Fe 3+ ions, respectively). We consider 
all 15 possible superexchange (SE) interactions and all 19 super-superexchange (SSE) in- 
teractions with the 0...0 distance less than 3.2 A. The intra- and inter-sheet interactions 
within each Fe204 layer as well as the SSE interactions between adjacent Fe204 layers are 
taken into account. To evaluate these 34 spin exchange parameters reliably, we considered 
111 different ordered spin states leading to 110 energy differences. The 34 spin exchange 
parameters were determined by performing a linear least-square fitting analysis. The SSE 
interactions are generally much weaker than the SE interactions with the magnitude of all 
SSE interactions less than 1.4 meV. The calculated SE parameters are reported in Table [B 
All intra-sheet SE interactions are AFM, and the strongest interactions (~ 7.3 meV) occurs 
between the 3B Fe 3+ ions because of the large energy gain of the AFM configuration and al- 
most zero FM coupling. The inter-sheet SE interactions are weaker than the the intra-sheet 
SE interactions, and are mostly AFM. 

With the calculated spin exchange parameters, one can identify the spin ground state of 
the CO state. The Metropolis Monte Carlo simulation of the Ising model is performed to 
search for the ground state. Simulations with supercells of several different sizes show that 
the spin ground state has the magnetic structure shown in Fig. |2](a), which has the same cell 
as the \/3 x v3 CO structure. In this state, all Fe 2+ ions contribute to the majority spin, 
and the Fe 3+ ions are antiferromagnetically coupled to the Fe 2+ ions in the type A T-sheet. 
In the honeycomb lattice of the type B T-sheet, the Fe 3+ spins are antiferromagnetically 
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coupled. Thus, the spin ground state is ferrimagnetic, as experimentally observed This 
2:1 ferrimagnetic structure is the same as the magnetic structure proposed by Siratori et al. 
[3], and differs from the structure proposed by Naka et al. 

The observed ferrimagnetic ordering can be readily explained in terms of the calculated 
exchange parameters. In the honeycomb network of the type B T-sheet, the nearest-neighbor 
(NN) 3B ions are antiferromagnetically coupled since their SE interaction is strongly AFM. 
In the type A T-sheet, the SE interactions between the 2A ions are AFM, and so are those 
between the 2A and 3A ions, which leads to spin frustration. As a consequence, two possible 
spin arrangements compete with each other in the type A T-sheet; the first is the state in 
which the coupling between the NN 2 A ions are AFM with the spin direction of the 3 A 
ion undetermined, and the second is the state in which all 2A ions are antiferromagnetically 
coupled to the 3A ions. The energies of these two states (considering only the SE interaction) 
are E 1 = -4(J 2 ai,2A2 + Jzaim*) per 3A ion, and E 2 = -10( Jza\,ia\ + <hA\,2Ai + J 3Ai,2As) + 
4(</2Ai,2A2 + J2Ai,2A4) per 3A ion, respectively. Due to the relatively strong AFM interactions 
between the 3A and 2A ions (See Table [T|) and the large spin of the 3A ions, the second 
state has a lower energy, i.e., E 2 < E\. Without loss of generality, we can assume the 2A 
(3A) ions constitute the majority (minority) spin in the second state. Now, we examine 
the spin orientation of the Fe 2+ ions in the type B T-sheet. The intra-sheet interactions 
of the 2B ion with 3B ions vanish due to the AFM ordering of the 3B ions. As for the 
inter-sheet interactions involving the 2B ions, the dominant one is the AFM interaction of 
the 2B ion with the 3A ion (J3A1-2B1 in Tabled]). Consequently, we obtain the ferrimangetic 
ground state shown in Fig. 12(a), in which the spin of the 2B ion contributes to the majority 
spin of the Fe 2 4 layer. For the stability of the ferrimangetic ground state, the inter-sheet 
interaction is essential. This was neglected in the model Hamiltonian study of Naka et al. 
jl8| . The ferrimangetic state is not due to the FM interactions between NNN Fe ions of the 
T-sheet because they must be vanishingly weak and mostly AFM. 

The electronic structure of the ferrimangetic state calculated for the v3 x v3 CO struc- 
ture of LuFe204 is shown in Fig. [3j Also shown is the electronic structure calculated for 
the FM state. Both states are semiconducting, and the highest occupied (HO) and the 
lowest unoccupied (LU) levels of both states come from the spin-up Fe 2+ and Fe 3+ ions, 
respectively [5[. In addition, the band dispersion from T to A is rather small, indicating a 
very weak interlayer interaction. However, there are some important differences. First, the 



5 



ferrimangetic state has a larger band gap (1.68 eV) than does the FM state (0.77 eV). This is 
consistent with the stability of the ferrimangetic state. Second, the FM state has an indirect 
band gap with the HO and LU levels located at K and T, respectively. In the ferrimangetic 
state, however, the LU level has the highest energy at T and the band dispersions of the HO 
and LU levels are almost flat from M to K. This difference comes from the orbital interaction 
between the spin down (d x 2_ y 2,d xy ) levels of the spin up Fe 3+ and Fe 2+ ions. 

To probe the presence of spin-charge coupling in LuFe204, it is necessary to consider the 
spin orderingin a CO state other than the a/3 x a/3 CO state. The previous electrostatic 
calculations flQ showed that the chain CO, in which one- dimensional (ID) chains of Fe 2+ 
ions alternate with ID chains of Fe 3+ ions in each T-sheet [Fig. 12(b)], is only slightly less 
stable than the a/3 x a/3 CO, and has no FE polarization. We extract exchange parameters 
by mapping analysis as described above. It is found that the intra-sheet SE between the 
Fe 3+ ions is the strongest ( J = 6.7 meV) as in the a/3 x a/3 CO case. All intra-sheet SE's are 
AFM with J(Fe 3+ -Fe 3+ ) > J(Fe 2+ -Fe 3+ ) > J(Fe 2+ -Fe 2+ ). The inter-sheet SE between the 
Fe 3+ ions is very weak (| J\ < 0.3 meV), and that between the Fe 2+ and Fe 3+ ions is FM with 
J = —1.4 meV. Interestingly, the inter-sheet SE between the Fe 2+ ions is rather strongly 
AFM (J = 6.3 meV). Monte Carlo simulations using these spin exchange parameters indicate 
that the spin state shown in Fig. [2(b) is the spin ground state. In this spin ordering, the 
spins within each chain of Fe 2+ ions or Fe 3+ ions are antiferromagnetically coupled. The 
NN chains of Fe 2+ ions belonging to different T-sheets are coupled antiferromagnetically, 
whereas the corresponding chains of Fe 3+ are almost decoupled. 

The above results show that the spin ordering of the chain CO state is dramatically 
different from that of the v3 x a/3 CO state. The most important difference is that the 
total spin moments are 2.33 /zb/FU for the a/3 x a/3 CO, but z/#/FU for the chain CO. 
This evidences a strong spin-charge coupling in LuFe204. The external magnetic field will 
have different effects on the two CO states due to the the Zeeman effect. It is expected that 
the magnetic field will further stabilize the ferrimagnetic a/3 x a/3 CO state. Consequently, 
an external magnetic field will reduce the extent of charge fluctuation and hence decrease 
the dielectric constant. This supports our explanation for the giant magnetocapacitance 
effect of LuFe 2 4 at room temperature [5( . 



Without considering the inter-sheet interactions, Naka et al. 



18| suggested that the 



degeneracy of the spin ground state of the a/3 x a/3 CO state is of the order 0(2 JV / 3 )( N is 
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the number of the spin sites), which is much larger than the spin degeneracy [0(2^)] of 
the chain CO state. Thus, they proposed that spin frustration induces reinforcement of the 
polar x ^3 CO by a gain of spin entropy. However, our calculations show that there 
are substantial inter-sheet spin exchange interactions between the 2B1 and 3A1 ions, which 
would remove the macroscopic degeneracy of the spin ground state of the y/3 x y/3 CO state. 
The macroscopic degeneracy still persists for the chain CO state. Thus, our work provides 
a picture opposite to what Naka et al. proposed. Furthermore, we find that the y/3 x y/3 
CO state is more favorable for the spin ordering than is the chain CO state; with respect to 
the paramagnetic state, the spin ground state is lower in energy by —78 meV/FU for the 
y/3 x y/3 CO, but by -57 meV/FU for the chain CO. The model of Naka et al. [18] predicts 
that the polar y/3 x y/3 CO state is destabilized and the electric polarization is reduced by 
the magnetic field, since it will lift the macroscopic spin degeneracy. In contrast, our work 
predicts that the magnetic field stabilizes the ferrimagnetic y/3 x y/3 CO state due to the 
Zeeman effect, and provides an explanation for why the electric polarization increases when 
the temperature is lowered below the Neel temperature jsj], because the charge fluctuation 
has an onset well below Tqo 

In summary, our first principles results explain the experimentally observed Ising ferri- 
magnetism, and manifest the spin-charge coupling and magnetoelectric effect in LuFe204. 

Work at NREL was supported by the U.S. Department of Energy, under Contract No. 
DE-AC36-08GO28308, and work at NCSU by the U. S. Department of Energy, under Grant 
DE-FG02-86ER45259. 
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TABLE I: Calculated superexchange parameters (in meV) in the v3 x ^/?> CO state of LuFe204 
(For the spin sites of the 2A, 3A, 2B and 3B ions, see Fig. Q] ) 
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FIG. 1: (Color online) Schematic representation of the v3x v3 CO structure. Large, medium, and 
small circles represent the Fe 2+ , Fe 3+ , and 2 ~ ions, respectively. The type A (type B) T-sheet has 
the honeycomb network of Fe 2+ (Fe 3+ ) ions with a Fe 3+ (Fe 2+ ) ion at the center of each hexagon. 
2A and 3A (2B and 3B) refer to the Fe 2+ and Fe 3+ ions of the type A (type B) T-sheet, respetively. 
The region enclosed by dashed lines indicates the unit cell of the CO structure. There is a mirror 
plane of symmetry, which is parallel to the c axis and crosses the 3A1 and 2B1 sites. The inset 
shows an isolated FeOs trigonal bipyramid. 



FIG. 2: (Color online) Schematic representations of (a) the spin ground state of the x \/3 
CO structure and (b) one of the macroscopic spin ground states of the chain CO structure. The 
arrows denote the spin directions. The region enclosed by the dashed lines on the bottom T-sheet 
indicates the magnetic unit cell of the spin structure. 
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FIG. 3: (Color online) Band structures calculated for (a) the FM state and (b) the ferrimagnetic 
state of the \/3 x \/3 CO structure of LuFe204. The solid and dashed lines represent the up-spin 
and down-spin bands, respectively. The \/3 x \/3 x 1 hexagonal cell is used in the calculations. 
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